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Abstract 

Metacommunity theory is considered a promising approach for ex- 
plaining species diversity and food web complexity. Recently Pillai et 
al. proposed a simple modeling framework for the dynamics of food 
webs at the metacommunity level. Here, we employ this framework to 
compute general conditions for the persistence of complex food webs 
in metacommunities. The persistence conditions found depend on the 
connectivity of the resource patches and the structure of the assembled 
food web, thus linking the underlying spatial patch-network and the 
species interaction network. We find that the persistence of omnivores 
is more likely when it is feeding on (a) prey on low trophic levels, and 
(b) prey on similar trophic levels. 
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1 Introduction 



A central aim of ecology is to understand the emergence and maintenance of 
the enormous diversity of ecological species. Both in aquatic and terrestrial 
environments the number of niches created by primary abiotic factors is 
relatively limited. Most ecological communities thus exhibit a self-sustaining 
diversity where species present in the system open up niches for others. 

Perhaps the simplest example for biotically created niches is found in 
predator-prey interactions. By persisting in a given system a species can 
open up a niche for a predator feeding on that species. This mechanism can 
in principle enable the coexistence of a large numbers of species in complex 
food webs. 

Recent investigations have identified several factors that contribute to the 
stability of large food webs and thus promote diversity (McCann, Hastings 
& Huxel 1998, Brose, Williams & Martinez 2006, Vandermeer 2006, Gross, 
Rudolf, Levin & Dieckmann 2009). However, growing evidence suggests that, 
at least in some systems, the food web emerges only on a regional scale, 
whereas simple food chains are observed if specific locations (patches) are 
considered in isolation (see Fig. [I]) (Pillai, Gonzalez & Loreau 2011). This 
points to a need for a met a- community perspective, in which one explicitly 
accounts for the dispersal of species across a network of patches. In par- 
ticular, one can ask if and how different communities can be sustained in a 
network of similar patches. 

The dynamics of dispersal in a system of multiple patches was first 
discussed in the context of metapopulation theory, introduced by Levins 
(Levins 1969, Hanski 1998). The simplest metapopulation model considers 
only a single species that spreads over a system of identical patches, such 
that each patch is either occupied by the species or vacant. Extending this 
framework to allow for multiple populations per patch leads to metacommu- 
nities. 

In the present paper we focus on a class of metacommunity models that 
has recently been proposed by Pillai et al. (Pillai, Loreau & Gonzalez 2010, 
Pillai et al. 2011, Pillai, Gonzalez & Loreau 2012). These models are rela- 
tively simple and therefore conducive to a detailed analysis. At their center 
stands the assumption that within each patch competitive exclusion pre- 
cludes the formation of complex food webs. Thus any single patch can only 
harbor a simple food chain. However, as different food chains can be realized 
in different patches a complex food web emerges at the landscape level (see 
Fig. §. 

Pillai et al. discuss the conditions for specific examples, when complex 
food webs can persist indefinitely in the metacommunity. Here we continue 
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Figure 1: Schematic illustration of the emergence of complex food webs. 
Spatially coupled local patches sustaining food chains give rise to complex 
food webs on the landscape scale. Figure reprinted from (Pillai et al. 2011) . 



this work and show that the computation of persistence conditions follows a 
general pattern. Exploiting this pattern, we present a general approach al- 
lowing the direct computation of conditions for the persistence of an arbitrary 
given food web in the metacommunity. 

We start in Sec. [2] with a detailed description of the modeling framework. 
In Sec. [3| we compute conditions for food chain persistence, introducing a 
concept which is then used to compute persistence conditions for more com- 
plex food webs in the subsequent sections. Based on Sec.|4| where we consider 
a simple example involving omnivory, we formulate a general formalism for 
the calculation of persistence ranges in Sec. |5j In Sees. [6] and [7] we use 
this formalism to investigate the dependence of the persistence range on the 
predator-prey interactions and on the geographical network. We conclude 
with a summary and discussions in Sec. [8j 

2 Modeling framework 

Following Pillai et al. (Pillai et al. 2010, Pillai et al. 2011) we consider a 
metacommunity consisting of s species that populate a network of discrete, 
interconnected patches. We do not account for the abundance in specific 
patches, but capture the presence (or absence) of a certain species in a par- 
ticular patch by a binary variable. Thus every individual patch is either 
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empty or occupied by a specific set of species. 

The system evolves due to local extinctions and colonization of adjacent 
patches. Local extinctions occur either spontaneously (e.g. due to external 
events or demographic stochasticity) or due to interactions with other species 
in the model. For every species % we denote the rate of spontaneous local 
extinctions as ej. When the extinction of a species occurs, all species directly 
or indirectly feeding upon this species also go locally extinct. For example, 
in a patch with three species where species 3 feeds on species 2 and species 2 
feeds on species 1, species 1 goes extinct with rate e±, species 2 goes extinct 
with rate t\ + e<i and species 3 goes extinct with rate e\ + e2 + e%. 

Colonization allows a species % that is established in a patch to establish 
itself in suitable patches with a constant rate q. A target patch is considered 
suitable if a) it can be reached from the source patch, b) no stronger com- 
petitor is already established in the target patch, and c) prey for the focal 
species is established in the target patch. 

Regarding a) we assume that the patches form a complex network, to 
which we refer as the geographical network. In this network every node rep- 
resents a patch, and a link between two nodes indicates that colonization is 
possible between the corresponding patches. We assume that the geograph- 
ical network is identical for all species, which is not necessarily true in all 
systems. Furthermore, we note that assuming constant colonization rates per 
link implies that dispersal is not strongly limited by the number of colonizers. 

Regarding competition, b) we assume that specialists are better adapted 
and therefore superior competitors than generalists. If there is an exploita- 
tive competition between a specialist and a generalist predator in the same 
patch and the generalist cannot feed on the specialist then the generalist is 
excluded. This is assumed to occur instantly in the model. Thus, a general- 
ist predator cannot colonize a patch where a competing specialist is already 
established. Furthermore, an established generalist population goes locally 
extinct if a competing specialist colonizes its patch. 

Regarding the availability of suitable prey, c), we assume that only certain 
primary species can colonize empty patches. All other species require suitable 
prey species in their respective patches. In any particular model system 
we thus choose a fixed set of suitable prey species for each species in the 
model. Together all of these potential predator-prey relationships also define 
a network. In this trophic network the nodes represent the species, whereas 
directed links represent potential predator-prey interactions. We note that 
every food chain observed in a particular patch is a subgraph of this trophic 
network (i.e. a subset of its nodes and links). Therefore, also the food web 
that is formed on the landscape scale, by superposition of the food chains 
found in patches, is a subgraph of the trophic network. The trophic network 
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is thus the maximal food web that can be observed on the landscape scale. 
However, the observed food webs may be missing some nodes, if some species 
go extinct on all patches, or links, if some predator-prey interactions are not 
realized in any patch. 

In the following sections we investigate the conditions for the coexistence 
of species on the landscape scale, depending on the structure of the trophic 
and geographical networks. 

3 Linear food chain 

We start by considering the case where the maximal trophic network is a 
linear tri-trophic food chain (see Fig. [2] a)). The main purpose of this simple 
example is to illustrate the calculation of coexistence ranges, which is ex- 
tended into a general rule for more complex cases in the subsequent sections. 
We note that a very similar system, although on a different network, was 
already analyzed in (Pillai et al. 2011). 

To gain a mathematical understanding of the metacommunity dynamics 
we formulate a mean-field model describing the density of patches in which a 
specific set of species is found. We use a notation where symbols of the form 
[i] indicate the proportion of patches where the local food chain has length 
i. Thus, [0], [1], [2] and [3] denote the proportion of patches that are empty, 
inhabited only by species 1, inhabited only by species 1 and species 2, and 
inhabited by all three species, respectively. Conservation of the number of 
patches implies [0] + [1] + [2] + [3] = 1. Additionally we use pi, p 2 and p 3 
to denote the proportion of patches where species 1, species 2 and species 
3 is present, irrespective of the presence of other species. Consequently, 
Pl = [1] + [2] + [3),p 2 = [2] + [3] and p 3 = [3]. 

The expected time evolution of the mean-field patch densities is then 
given by 

£[0] = -c 1 (k)[0]p 1 + e& 1 

|[1] = Cl (k)[0} Pl - c 2 (k)[l]p 2 + e 2 p 2 - e x [l] 

| [2] = c 2 (k)[l]p 2 - c 3 (k}[2}p 3 + e 3 p 3 - ( ei + e 2 )[2] 

|[3] = C3(fc)[2]p3-(e 1 + e 2 + e 3 )[3], 

where we introduced (k) to denote the mean number of neighboring patches 
per patch. 

In Eq. the first terms on the right-hand-sides describe the effect of 
colonization, whereas the second terms describe the effect of extinctions. To 
understand the specific functional form appearing in the equations let us 
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consider the first term of the first equation, —ci(k)[0]pi, which describes 
the loss of empty patches due to colonization. We lose empty patches due to 
colonization from adjacent patches in which species 1 is established. A typical 
patch has (k) adjacent patches and species 1 is established in any of these 
with probability p±. Multiplying the expected number of adjacent patches 
where species 1 is established, pi(k), with the colonization rate c\ yields the 
total expected colonization rate for a given empty patch. Multiplying further 
with the proportion of empty patches, [0], yields the total loss rate Ci(k)[0}pi. 

We note that in the equations for empty patches ([0]) and tri-trophic 
patches ([3]) the effect of colonization is purely negative and purely positive, 
respectively. By contrast, for single-species ([1]) and two-species ([2]) patches 
colonization appears both, as a gain and as a loss term, as patches of these 
types can be both, created and destroyed by colonization. Similarly, the 
effect of extinctions is purely positive for [0]-patches and purely negative for 
[3]-patches, whereas for [1]- and [2]-patches extinction has both, positive and 
negative effects. 

In the equations we replaced the number of patches adjacent to a given 
patch, with the expectation value of adjacent patches (k), the so-called mean 
degree. Working with this mean degree instead of the full probability distri- 
bution of the numbers of adjacent patches (the degree distribution) is known 
as homogeneous approximation in network science. This approximation is 
known to yield good results for all reasonably narrow degree distributions, 
which should cover most ecologically relevant cases. However, we emphasize 
that a somewhat more complicated computation is necessary for instance for 
geographical networks with a central hub. 

By setting the right-hand sides of Eq. ([TJ to zero and solving the resulting 
system of equations, we find the expected equilibrium patch densities 
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where we used asterisks to indicate the stationary values that the system 
approaches after sufficiently long time. 

Generalizing these results to a class of systems where the maximal net- 
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work is a linear food chain of length Z, one finds 



e x H h e i+ i . 

— , < i < I 

c i+ i{k) 



w = i - 



i=0 



In the following we use similar reasoning to study more complex cases. To 
simplify the presentation we assume that the extinction rates for all species 
are equal, = e Vi, and further that the colonization rates for all species are 
equal, c, = c Vz. Additionally, we introduce the dimensionless parameters 
z = e/c and z = z/(k), which constitute differently normalized extinction 
rates. In this notation the equilibrium patch densities are 

= (i + 1)2, < i < I 

[1}* = 1 - + l)f , (2) 
p? = 1 - §i(z + 1)^, l<i<l. 

Let us now investigate under which conditions a food chain of length 
/ can persist in the metacommunity. To obtain an answer directly, we 
use an approach motivated by the epidemic literature (Kamgang & Sallet 
2008, Van den Driessche, Watmough et al. 2002, Diekmann, Heesterbeek & 
Roberts 2010, Bohme & Gross 2011) and compute the linear stability of the 
"Z-free equilibrium", i.e. the equilibrium where the species / is absent, which 
is analogous to considering the "disease-free equilibrium" for epidemics. In 
other words, we analyze the stability of a food chain of length / — 1 against 
arrival of species /. This reveals the threshold parameter values at which 
the / — 1-trophic food chain becomes unstable and species / can establish 
itself. As will become apparent below, this is also the threshold for persis- 
tence of species / because the present model cannot exhibit Allee-effects on 
the metapopulation level. 

From the evolution equation for the patch density [I], 

f t [l] = (c(k)[l - 1} - le)[l], 

it follows that persistence of species I is possible if 

{c(k)[l- 1]* -le)[l] > 0. (3) 

We note that [/ — 1]* denotes the equilibrium patch density of species I — 1 
in the "Z-free equilibrium", before arrival of species I. The condition Eq. ^ 
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shows that the prey, species I — 1, must reach a threshold density [I — 1]* = Iz 
to allow for persistence of species I, which confirms basic ecological intuition. 
Using Eq. ^ the persistence condition for I can be rewritten as 

* < (4) 

which can be read as a bound for the dimensionless extinction rate z. Thus, 
high extinction rates hinder the formation of long food chains while a dense 
geographical connectivity (high (k)) promotes it. 



4 Omnivory in a tri-trophic example 

In this section, we calculate the parameter range for z where food webs 
containing specialists and omnivores can persist at the metacommunity scale. 
These analyses generalize the results of (Pillai et al. 2011) by providing a 
general formalism for the computation of persistence ranges. 

As a first example, we consider the food web configuration shown in 
Fig. [2j Here, the maximal food web consists of a linear food chain and an 
additional omnivore top predator that can feed on both the producer and the 
primary consumer in the food chain. We investigate under which conditions 
the omnivore, x, can coexist with the specialist top-predator of the food 
chain, although the specialist will exclude the omnivore in any given patch. 

Even without mathematical reasoning it is intuitive that the omnivore 
faces two dangers. First, if the geographical network is very sparse then the 
prey species may be too sparse to allow for persistence of the omnivore. Sec- 
ond, if the geographical network is too dense then the specialist top-predator 
will establish, which threatens the omnivore with competitive exclusion. 

In the following we mathematically determine the threshold connectivities 
of the geographical network by employing the same line of reasoning as for 
the linear chain. We start by writing the governing equations for patches in 
which the omnivore x is established as 

IN = c(k)([l)p x - [lx)p 2 ) + e([2x) - 2[lx)) 
f t \2x] = c(k)([2] Px + [lx]p2 - [2x]p 3 ) - 3e[2x], 

where we used symbols of the form [ix] to denote patches in which a specialist 
food chain of length i and the omnivore x is established. Note that omnivore x 
does not suffer competitive exclusion in [lx]-patches, because species 2 is not 
only a competitor for x, but also one of its prey species. Thus, colonization of 
species 2 turns a [lx] -patch into a [2x]-patch. This contribution is captured 
in the term c(k)[lx]p2- 
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Figure 2: Example for a maximal food web emerging in the metacommunity. 
Specialist species are depicted by black nodes, omnivores by gray nodes. 
Specialists are numbered according to the trophic level they belong to. Ar- 
rows indicate feeding links. Shown are different food web that arise on the 
landscape-scale in the parameter rangesindicated. c) The omnivore x is able 
to persist in a system of two specialists for z < {k)/3, coexists with the 
specialist species 3 for z < (k)/6 and becomes extinct due to competitive 
exclusion for z < (k)/(6 + 3a/2). These parameter ranges are not drawn to 
scale. See text for derivation. 



To allow for a more concise treatment we can rewrite Eqs. ^ in the 
matrix form 




/ c(k)([l]-p 2 )-2e c(k)[l]+e \ /[lx] 
c(k)([2]+p 2 ) c(k}([2]-p 3 )-3e) \[2x] 



(6) 



s 



describing the time evolution of the vector p x = ([lx], [2a;]). 

Now, we analyze the response of the x-free equilibrium to small pertur- 
bations, i.e. arrival of a; in a small number of patches. Persistence of the 
omnivore x is possible when the x-free equilibrium is unstable. From nonlin- 
ear dynamics we know that this is the case when the largest eigenvalue A of 
the matrix S, in the a>free equilibrium, is positive, which we write as 

A(S*) > 0. (7) 

Starting from this condition we now extract an explicit condition in terms of 
the connectivity of the geographical network (k) and the extinction rate z. 

In the evaluation of the eigenvalue inequality, we have to take into account 
that S* depends on the equilibrium patch density of species 3. Because it 
thus makes a difference whether species 3 is present or absent, we have to 
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analyze the stability of the x-free equilibrium separately for these two cases 
(see Fig. [2]b) and c)). To distinguish between the two cases, we use S* to 
denote S* when species 3 is present (0 < z < (k)/6) and S* to denote S* 
when species 3 is absent ((h)/ 6 < z < (k)/3). The persistence range for 
the omnivore x is then obtained from a combination of the two conditions 
A(S*) > and A(S*) > . 

Let us first consider the case when species 3 is present. Inserting the 
corresponding equilibrium expressions from Eqs. @ into Eq. g, we can 
write 



si s 2 



The matrix S* is thus split into two matrices, s x and s 2 , which capture the net 
balance of x-patches due to colonization and local extinction, respectively. 

From ^ it is clear that the x-free equilibrium is stable at z = 0, because 
all eigenvalues of S* =0 = s x are negative. When z is increased a threshold z' 
is eventually reached where one eigenvalue becomes positive and the omni- 
vore can no longer persist. Mathematically we can detect the transition by 
monitoring the determinant det(S*), which becomes negative at z' . In order 
to calculate the threshold, we write 

det(S*) = det(Si(Jfe) + s 2 z) = det(-Bi) det(Mz - (Jfe)), (9) 

where M = — s 2 sj L 1 , and s^ 1 is the inverse matrix of §i. From Eq. ^ it 
follows that det(S*) < if either a) det(— §i) < and all (real) eigenvalues of 
M smaller than (k)/z, or b) det(-Si) > and A(M) > (k)/z, where A(M) is 
the largest (real) eigenvalue of M. As for the matrix given in Eq. (|8]) we have 
det(— Si) > 0, the omnivore persists in the system for z > z' = (fe)/A(M). 

Now we have an explicit expression for the threshold z' . For the matrix 
S* in Eq. @, 

'3 3\ fl 0\ _ (6 3 N 
.0 6 J "VI I ~ U 6 



M 



and therefore, the persistence condition for p 3 ^ becomes 

l> 5^5- (10) 

This is a lower bound for the parameter z, which limits the persistence of x 
due to competition with the specialist 3. 
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Now, we consider the case p% = 0, shown in Fig. [2] c). Inserting the 
corresponding expressions from Eqs. ^ in Eq. ^ yields 



Si s 2 



In this case, we observe that the stability of the x-free equilibrium is lost when 
approaching some threshold value z' from above. This can be understood 
intuitively, because the x-free equilibrium is certainly stable when there is 
no prey species 2, i.e. for z > (k)/3. It becomes unstable if one eigenvalue of 
S* becomes positive at a value z', and accordingly det(S*) becomes negative. 
An analogous argumentation as before in Eq. ^ leads to the conclusion that 
A(S*) > is fulfilled when z < z' = (k)/X(M), where M = -s 2 s^ and we 



used that det(— §i) < for the matrix given in Eq. (11). 



For the matrix S* in (11) the corresponding condition yields 



(k) , v 

* < V- (12) 

This constitutes an upper bound for z, below which omnivore x persists in 
a system of two specialists. Note that we obtain the same upper bound for 
the persistence range of a specialist species 2, according to Eq. Q. Thus, as 
soon as a chain of two specialists exists in the metacommunity, an omnivore 
feeding upon both of them can establish. Below we show that this holds for 
any omnivore feeding upon two prey species on subsequent trophic levels. 



Taking conditions Eq. (10) and Eq. (12) together, we conclude that for 



<*> <:< ( 4 (13) 



6 + 3^ 3 

omnivore x is able to coexist with a chain of specialists at the metacommunity 
level. 

A summary of the results for this example system are shown in Fig. [2} 
We found that for z < (k)/3 omnivore x can establish itself in the meta- 
community. In contrast to the persistence range of species 2, the persistence 
range of omnivore x is bounded from below, due to competition with spe- 
cialist 3, which can persist in the system for z < {k)/6. Omnivore x can 
coexist with the superior competitor 3 for z > (k)/(G + 3\/2) (b), while for 
z < (k)/(6 + 3 a/2) the omnivore becomes extinct and the observed food web 
is a linear chain of three specialists (a). 

The results obtained in this section reproduce those already obtained by 
(Pillai et al. 2011). Whereas (Pillai et al. 2011) obtained these results from 
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1. Construct matrix S from evolution equations for p* x . 

2. Insert the corresponding equilibrium densities from Eqs. ^ in order 
to obtain S* and S*. 

3. Write S* in the form Si(fc) + s 2 z. Do the same for S*. 

4. Determine M = — S2Sj L 1 and A(M). Do the same for M. 

5. The persistence condition is (&)/A(M) < z < (k)/\(M). 



Figure 3: Recipe for the calculation of persistence thresholds. 

the solution of the full underlying equation system, we compute them directly 
from a stability analysis. The advantage of this approach is that it can be 
easily applied to more complicated situations and lends itself to automation 
in computer algebra systems, which allows to compute coexistence ranges for 
a large variety of different food webs. 

5 General method for persistence thresholds 

We now exploit the power of the method proposed in the previous section to 
obtain general results on the persistence of omnivores. Consider an omnivore 
feeding on different prey species, where the trophic level of the highest prey 
species is i — 1. In this case the corresponding S*-matrix is given by 

2(k) 



S* = §i(fc) + s 2 z } for 0<z< 



S* = h(k)+s 2 z, for < z < 2{k) 



i(i + 1) ~ ~ 1) 

The different ranges appear in this equation because of the extinction of 
species i and extinction of species i — 1, respectively. From the general 
expression for S* we can proceed analogously to the previously discussed 
specific example. Importantly, certain properties of the specific example sys- 
tem that we exploited above remain valid in the general case. For example, 
S* has only negative eigenvalues at z = 0, because Si is always lower trian- 
gular and has negative diagonal entries. Therefore det(— Si) > 0, and due 
to Eq. @ there is a lower bound for the persistence range, which is deter- 
mined by A(M). Conversely, as det(— §i) < 0, there is an upper bound of 
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the persistence range, if A(M) = (k)/z in the relevant range. Otherwise, the 
persistence range is naturally bounded at 2(k)/(i(i — 1)), where the highest 
prey species of the omnivore (species i — 1) goes extinct. 
Summarizing the considerations above, we find 

A(S*) > in the range — t~^~t <*< (14) 
V ; A(M) A(M) V ; 

This relationship enables the direct computation of the persistence range for 
a given omnivore, and thus allows for the analysis of a broad class of food 
webs in metacommunities. An overview of this procedure is given in Fig. [3} 
An additional example for the application of this procedure to a specific food 
web is given in the Appendix. 



6 General results on omnivory 

Let us now use the method proposed above to gain ecological insights into the 
persistence of omnivores in food chains. For simplicity, we restrict our anal- 
ysis to omnivores feeding upon two prey species, but note that the proposed 
method can be likewise applied to omnivores feeding on more prey species. 
Unless stated otherwise we denote the trophic levels of the omnivores' prey 
by i — 1 and j — 1, where i > j. Correspondingly, the specialists competing 
with the respective omnivore occupy trophic levels % and j. 

We distinguish between the following two cases (see Fig. I) The 
omnivore feeds upon prey at subsequent trophic levels (i — j + 1). II) The 
omnivore feeds upon prey at non-subsequent trophic levels (i > j + 1). The 
S*-matrices governing the evolution equations for [(i — l)x}- and [(j — l)x}- 
patches in the x-free equilibrium are then given by 

f(k)(\j-l]*-p*)-jz (k)\j-l]* + z \ 

Si* = (15) 

V (k)([i-l]*+ P U) (*)([« -1]* -#)-«/ 

and 

f(k)([j-l]*-p*)-jz (k)[j-l]* + z \ 

Sn = , (16) 

V (k)[i-i\* (*>([< -l]* -#)-<*/ 

where we introduced the subscripts I and II in order to distinguish the S*- 
matrices in the different cases. 

In order to calculate the coexistence range for a general omnivore, we 
have to separately consider the situation where species i is absent (p* = 
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Figure 4: General food web configuration consisting of a chain of specialists 
and an omnivore x. In I) the omnivore x feeds on two subsequent trophic 
levels, i.e. i = whereas in II) the omnivore x feeds on two non-subsequent 
trophic levels, i.e. % > j + 1. Accordingly, different coexistence ranges are 
obtained for both CclSGS, clS described in the text. 



0, [i — 1]* = p*_ 1 = 1 — — l)z)/2), and the situation where species i is 
present in the system (p* = 1 — + l)z)/2, [i — 1}* = iz). Then, combining 
both conditions, we expect to find a parameter range for z which is bounded 
from below and from above. 



Following the procedure in Fig. |3|and using Eq. (|14|), we find 
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and the corresponding coexistence ranges 

(*> . . W 



I) 
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A (Mn) 



< z < 



< z < 



A(M : ) 

(k) 
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(17) 
(18) 



Equations ( 17) and ( 18 ) reveal the impact of the connectivity of the patch- 
network, showing that dense topologies promote the persistence of omnivores 
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Figure 5: Dependence of the coexistence range for an omnivore on the trophic 
levels of its prey. The grey regions show the coexistence ranges for different 
values of % and varying j, corresponding to an omnivore feeding on species 
i — 1 and j — 1 (see Eqs. (17), (18)). The coexistence range becomes smaller 
with increasing % and for fixed % with decreasing j . This shows that omnivores 
are more likely to coexist if they feed on adjacent tropic levels, low in the 
food chain. Dashed lines indicate the respective persistence thresholds of 
species i. The dashed region marks the parameter range where the maximal 
number of omnivores coexists with the specialist chain. 
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in the metacommunity. We can further deduce statements concerning the 
dependence of the coexistence range on the trophic levels of the prey species: 



1. The coexistence range for omnivore x becomes smaller with increasing 
trophic level of the prey species. 

2. The coexistence range for omnivore x becomes smaller with increasing 
distance between the trophic levels of its prey species. 

These characteristics of the coexistence range are shown in Fig. [5j Note that 
if we are only interested in the range where the omnivore coexists with the 



competitor i, the upper thresholds in Eq. (17) and Eq. (18) become the per- 
sistence threshold for species i, i.e. z < 2(k)/(i(i + l)) (dashed line in Fig. [5J. 
The range where possibly more than one omnivore coexists with a chain of 
specialists is given by the intersection of the respective individual coexistence 
ranges. The dashed areas in Fig. [5] mark the parameter ranges where all co- 
existence ranges for the same value of i overlap. Thus the maximal number 
of coexisting omnivores in these parameter ranges yields n™ ax = % — 1. 

Let us emphasize that the coexistence range of an omnivore only depends 
on the trophic levels of its prey species and not on the total length of the 
specialist chain. This is a consequence of the following relationship: 

[i - 1] + [(i - l)s] = (f 5 (19) 

[l-p{i-l)z; Pi = 0. 

In the x-free equilibrium, i.e. for [(j — l)x] = [(i — l)x] = 0, the equilibrium 
densities given in Eq. ^ are recovered. Relation (19) can be interpreted as 



"biomass-conservation" , implying that the species densities in the chain 
are conserved, no matter how many predators feed upon species % — 1 (and 
upon prey at trophic levels below % — 1). 



The conservation law in Eq. (19) has another practical implication: it 
allows for the direct computation of equilibrium patch densities p x for om- 
nivores. In the Appendix we derive that the total density of patches where 
omnivore x is established is given by 

P* x =[(j-l)x)* + [(i-l)x)* = \(S*), (20) 

where S* is the corresponding transition matrix in the x-free equilibrium 



(matrices given in Eq. (15), Eq. (16)). Note that S* is altered at the point 
where species i enters the system. Thus, there is no unique matrix S* de- 
scribing p* x in the whole coexistence range. This is apparent in Fig. [6] at 
the parameter value z = 2/(i(i + 1)). So, from the point where invasion of 
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Figure 6: Equilibrium densities of three different omnivores from numeri- 
cal simulations (symbols) and theoretical predictions (lines). Dashed lines 
indicate parameter values where the respective specialist competitor goes 
extinct in the metacommunity. The equilibrium patch densities observed in 
simulations display a close to linear behavior and are well described by the 
analytical formula in Eq. (20). The parameters for the simulation are the 



same as in Fig. 10 



17 



A) : 


= B ) | 














1 


hi 


/l 2 /ll 





Figure 7: General food web configuration consisting of a chain of specialists 
and a generalist g. The generalist feeds on two prey species in two indepen- 
dent chains arising from two types of habitat. In A), the trophic level of both 
prey species is the same (i = j), in B) the trophic levels differ at least by 
one (i > j). The corresponding coexistence ranges for both configurations 
are derived in the text. 



an omnivore becomes possible, at the upper bound of the coexistence range, 
its patch density increases up to the point where the specialist competitor 
becomes able to persist in the system. At this point (dashed lines in Fig. [6]), 
the patch density of the omnivore starts to decrease due to competition with 
the specialist i. Both, increase and decrease proceed almost linearly, as it 
can be seen in Fig. [6j For i > j + 1 and p* ^ the linearity is exact because 
S* n = -1 + M n z. 

The analytical expression for the equilibrium patch density of omnivores 
is useful for the computation of persistence conditions of predators that are 
capable of feeding on other omnivores. Then, the S-matrix contains patch 
densities of omnivores, in addition to patch densities of specialists. 

7 Competition between generalists and spe- 
cialists 

So far, we considered all empty patches to be equal. In this section we 
consider a system which comprises two different types of habitat h± and hi, 
which we assume to occur with equal probablity. Primary producers are 
assumed to specialize on a specific type of habitat. Thus, each habitat type 
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sustains a distinct chain of specialists (see Fig. |7j). Now, we can ask for the 
conditions under which a generalist, feeding on prey from both chains, is able 
to persist in the metacommunity. 

For simplicity, we focus again on generalists feeding on exactly two prey 
species. In analogy to the previous analysis we denote the trophic levels of 
the prey species of the generalist g with i — 1 and j — 1, where j < i. The 
corresponding matrix S*, describing the evolution of [(j — l)g]- and [(i — l)g]- 
patches in the g-free equilibrium is then given by 

f(k)(\j-l)*-p*)-jz (k)[j-l)* \ 

S* = • (21) 

V (k)[i-l\* (k)([i-l]*-p*)-izj 

In contrast to above, [(j — l)g]-patches cannot become [(z — l)g]-patches or 
vice versa, as both prey species are from different, independent chains. 

For the calculation of the coexistence range, we distinguish the two cases 
(see Fig. [7]): A) The two prey species are at the same trophic level. B) The 
two prey species are at different trophic levels. This is necessary because 
generalists feeding upon two prey species at the same level, say i — 1, expe- 
rience no competition with specialists if species i is not present. In contrast, 
a generalist feeding upon one prey at level i — 1 and a second one at some 
lower level j '• — 1, even in the absence of specialist i, experiences competition 
from the specialist j. 

Proceeding according to Fig. |3j we find the coexistence ranges 

(k) (k) 

l(l + 6) l A 
(k) (k) 

B) -^t < z < ' 



A(M) A(M) 

In case A, we wrote the largest eigenvalues of the corresponding M-matrices 
explicitly, as they assume a very simple form. In case B, the threshold values 
are given through the following matrices, 



2i i(i + l)J ' \ 2i i(i + 1 

Here we used that the equilibrium patch densities for a two-habitat system 
are given byp* = 1/2 — i(i + l)z/2 and [j] = (j + l)z, for j < i. The change in 
Pi compared to a single-habitat system (see Eq. ([2])) is due to normalization. 

Analyzing the dependence of the coexistence ranges for generalists with 
respect to the trophic levels of their prey, we find qualitatively the same char- 
acteristics as for omnivores: The coexistence range increases with decreasing 
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trophic levels and with decreasing distance between the levels. One exception 
is the case where both prey levels are equal (i = j). Here, the coexistence 
range is smaller than in the asymmetric case i — j + 1, though still larger 
than for i > j + 1 (see Fig. [sj) . Thus, as above, the maximal coexistence 
range for generalists is found when feeding on adjacent trophic levels. 

In order to relate the coexistence range to the omnivores discussed above, 
we have to determine the coexistence range for omnivores in the presence of 
a second independent chain. One finds that the coexistence range for omni- 



vores, given in its general form in Eqs. (17), Eq. (18), aquires an additional 



factor of 1/2 in every threshold (see Appendix). Then, a comparison shows 
that the coexistence range for omnivores is always (slightly) larger than the 
corresponding one for generalists considered in this section (see Fig. [8]). 
Using the conservation equation (19), one can obtain an analogous re- 



sult for the patch density of generalists, as the one provided in Eq. (20) for 



omnivores. 



8 Summary and outlook 

In this paper we studied extensions of the metacommunity model proposed in 
(Pillai et al. 2011). We developed a mathematical method for the computa- 
tion of persistence and coexistence conditions for different species, based on a 
successive linear stability analysis of the system in equilibrium. This method 
allowed convenient computation of the range of connectivities of the under- 
lying spatial network, for which competing species can coexists in different 
scenarios. 

Among other results we used the proposed method to show that that the 
coexistence range is increased when a) predators feed on low trophic levels 
and b) when the distance between the trophic levels of their prey is small. 
This was found to be true both for predators feeding on species from the 
same chain of specialists and predators feeding on prey from different chains 
of specialists. 

For simplicity, we considered only species which are capable of feeding at 
most on two prey species. We emphasize, the mathematical method proposed 
here can be directly extended to predators feeding on a larger number of prey 
species. 

While the underlying model makes strong simplifying assumptions, the 
results obtained here and in the previous works (Pillai et al. 2010, Pillai 
et al. 2011, Pillai et al. 2012) seem plausible and thus probably merit fur- 
ther investigation in the future. A clear advantage of the simple model is 
that it is analytically tractable and thus offers an alternative to numerical 
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Figure 8: Comparison between the coexistence ranges for omnivores and 
generalists. For a system of two types of habitat and 10 trophic levels (20 
specialists in two branches), we plot the coexistence range for an omnivore 
x, feeding upon prey species from exclusively one branch and the coexistence 
range for a generalist g, feeding upon two prey species from both branches 
(inset). Throughout the whole j-range, the coexistence range for the omni- 
vore is larger than the coexistence range for the generalist. 
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investigations. This advantage is reinforced by the method proposed here, 
which allows direct computation of persistence conditions even in complex 
scenarios. 

A possible extension not discussed here is to consider directed geograph- 
ical networks. For instance for dispersal in river networks the directed flow 
between patches should not be neglected, which leads to a directed connec- 
tivity (Muneepeerakul, Weitz, Levin, Rinaldo & Rodriguez-Iturbe 2007). In 
this number of adjustments have to be made, for instance replacing 

the mean degree by the mean in-degree in the equations. Another promising 
direction for future research is to extend the present work to strongly hetero- 
geneous networks, such as scale-free networks, where the degree distribution 
cannot be approximated by the mean degree 

Finally, we note that the method proposed here may also be applicable 
in other systems. For example, the model discussed here is closely related 
to ongoing investigations on simultaneous spreading of several (competing) 
diseases in social networks (Karrer & Newman 2011, Newman 2005, Masuda 
& Konno 2006, Ahn, Jeong, Masuda & Noh 2006). 
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Appendix A: Systems with multiple types of 
habitat 



Starting from an empty patch network with equal distribution of h%- and 
/12-patches, specialists feeding on one or the other resource will arrange in 
two independent chains. The equilibrium patch-densities for empty habitat 
and / specialists feeding upon each other in each of the two chains are then 
given by 

W = [h 2 ]* = z 

= (i + l)z, \<i<l 

1 _ I 

2 2 



Following the same lines as in the case of one habitat type, the invasion 
condition for a specialist I becomes 



1(1 + 1)' 

Comparing this result to Eq. Q, we find that the invasion threshold is di- 
vided by two, due to the presence of two independent specialist chains. It is 
straight forward to verify that in general the existence of m different, equally 
distributed habitat types leads to the formation of m linear chains where the 
inequality 

ml(l + l)' 

determines the invasion threshold for specialist species at trophic level /. 
Consequently, the patch density pi is then given by p* = 1/m — i(i + 1)^/2. 
Analogously, it can be shown that the persistence thresholds for omnivores 
and generalists acquire an additional factor 1/m. 



Appendix B: Equilibrium patch density of om- 
nivores 

The evolution equation for the patch density p x of an omnivore x feeding on 
level j — I and i — 1 is given by 

f t p x = S x p x , (22) 
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where S x denotes the transition matrix for x-patches away from the x-free 
equilibrium andp x = ([(j — l)x], [(i— l)x]). For example, in case II (i > j + l), 
the corresponding transition matrix is given by 

'\j - i]-Pj-jz \j-i]+z 

[i - 1] [i-l)-pi- iz 



In contrast to the matrix in Eq. (16), describing the x-free equilibrium, the 
amount of suitable resource patches form omnivore x is reduced, as patches 
where the omnivore is already established are not available. However, with 
Eq. ( 19 ) we can relate the transition matrix S x to the corresponding matrix 
S* of the x-free case: 

c c _ i [0* - i)^] [U - 



Rewriting Eq. (22) reveals, 
S* - 



dtP x 



[(j-l)x\ [(j-l)x\ 

[(i-i) x \ [d-i)x\ 



Px)Px- 



So, at equilibrium (with x), the density of x-patches is given by 

p: = [(j-i)*r+[(i-i)*]* = A(s*). 



Appendix C: Example for omnivory in a 4- 
trophic chain 

Here, we consider an additional example to the one provided in the main 
text, in order to demostrate the procedure stated in Fig. [3] In this example, 
an omnivore feeds on two prey species which are not at subsequent trophic 
levels (see Fig. |9]). Consequently, for step 1) we have to use the S*-matrix 



given in (16) for % = 4 and j = 2, 

'(k)([l]*-p*)-2z (k)[l]* + z 

(k)[3}* <*>([3]*-pJ)-4*. 



S* 
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Figure 9: Example for a maximal food web emerging in the metacommu- 
nity. In contrast to Fig. [2j trophic levels of the omnivores' prey species 
differ by more than one. Symbols are the same as in Fig. [2j Shown are 
different food web configurations which arise in the given parameter ranges. 
The corresponding parameter ranges (not true to scale) are obtained by the 
mathematical formalism presented in the text. 



Following step 2) and step 3), we construct 



SI S 2 



~1 1 )<*>+( -6 -10 J ^' ^ = °' 



Sl s 2 



using the equilibrium patch densities in Eq. Then, in step 4), we deter- 
mine the matrices 



M 



3 3 \ (1 ON _ (3 3 

4 10 ) ' lo 1/ ~ U 10 



§ 2 -s" 1 



and 

M 



S2 

and calculate their largest eigenvalues: 




-3 
4 10 



A(M) = 13 + V97, A(M) = 5 + a/13. 
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Figure 10: Comparison of the predicted ranges with simulation results. Plot- 
ted are the equilibrium patch-densities (•), pi (o) and p*(x), obtained 
from network simulations. Eeach point corresponds to an average over 10 
independent simulation runs. We used iV = 10000, (k) = 20, c = 0.05 and 
varied z = e within the indicated range. Dashed lines correspond to the 
persistence thresholds predicted by our analytical calculations. Coexistence 
range for omnivore x: 2/(13+ v97) < z < l/(5 + \/l3), persistence threshold 
for specialist 3: z — 1/6 and persistence threshold for specialist 4: z = 1/10. 



According to step 5), the parameter range where an omnivore feeding on 
trophic levels 1 and 3 can coexist with a specialist chain is given by 



13 + V97 5 + V13 



Comparing this range to the one given in Eq. (13), one can observe that 
the whole range is shifted to smaller z-values, and that the total size of the 
coexistence range is smaller compared to the previous example. These two 
features resemble the general behavior shown in Fig. |5j 

A summary of the obtained results for this example system can be found 
in Fig. [9j In Fig. [10] we show simulation results for this system. It can be 
seen that the analytically estimated ranges correspond to disappearance or 
appearance of the respective species types, indicating that the predicted food 
web configurations actually are realized in the metacommunity. 
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